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Abstract 

Within the framework of the semiclassical approximation, we derive the Boltzmann equa- 
tion describing the dynamics of colorless plasmons in a hot QCD plasma. The probability 
of the plasmon-plasmon scattering at the leading order in the coupling constant is ob- 
tained. This probability is gauge-independent at least in the class of the covariant and 
temporal gauges. It is noted that the structure of the scattering kernel possesses impor- 
tant qualitative difference from the corresponding one in the Abelian plasma, in spite 
of the fact that we focused our study on the colorless soft excitations. It is shown that 
four-plasmon decay is suppressed by the power of g relative to the process of nonlinear 
scattering of plasmons by thermal particles at the soft momentum scale. It is stated that 
the former process becomes important in going to the ultrasoft region of the momentum 
scale. 
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1 Introduction 



For about three decades there has been an increasing interest in theoretical research into 
various dynamical properties of (ultra)relativistic many-particle systems. It is connected 
with manifold applications to various problems in astrophysical systems, modern cos- 
mology, in multiparton processes in experiment with high energy heavy ion collisions etc. 
The kinetic phenomenons, having a purely collective character, are one of the more impor- 
tant aspects of complicated dynamics of many-particle systems under extreme conditions. 
Here, the basic element in the description of transport phenomena is the derivation of 
the corresponding kinetic equations which would take into account (depending on the 
character of the problem being studied) the presence of mean fields in the system, two - 
(and more) body collisions, the possible renormalization effects, the effects of quantum 
fluctuations (stochastic), pair production, etc. Here, we restrict our consideration to a 
brief review of the results of the derivation of relativistic kinetic equations essentially 
based on two-body collisions in hot gauge theories. 

At present there are a few methods of construction of relativistic collision integrals. 
In particular, we mention the Zubarev's method of the non- equilibrium statistical operator 
[1] (use of this method on a relativistic systems of quark-gluon plasma (QGP) type can 
be found in Ref. [2]), and the method developed by Klimontovich [3] for an ordinary 
nonrelativistic plasma (the so-called second momentum or polarization approximation) 
and expanded to relativistic (semi) classical systems - in [4]. It should be stressed that 
the above-mentioned methods are particularly effective in the construction of collision 
integrals for relativistic (semi) classical systems, whose evaluation is described by so-called 
exact "microscopic" dynamical equations arising from the classical equations of motion. 
However, the extension of these methods to the relativistic quantum systems encounters 
some difficulties and therefore is uneffective. For the latter systems, mention may be 
made of the method based on the use of the reduction formulae of quantum field theory 
given by de Groot et al [5]. A more powerful and more convenient tool to derive the 
approximate relativistic kinetic equations from exact field Shwinger-Dyson equations, is 
the so-called closed-time-path (CTP) formalism [6]. Examples of its relativistic field- 
theoretical generalization can be found in [7, 8]. 

As is known, the cornerstone of derivation of the kinetic equations for a hot non- 
Abelian plasma is a fundamental separation of the momentum scale. The physical justi- 
fication for such a separation is the fact that the collective excitations which develop at 
a particular energy scale gT (or g 2 T, where T is the temperature, and g is the coupling 
constant), are well separated, when g -C 1, from the typical energies of plasma particles 
~ T [9]. Generally speaking one can define two types of kinetic equations: the equations 
for hard [] particles - hard quarks(g), antiquarks(g) and hard transverse gluons(g), and 

For the connection with our previous papers [10] here, we shall follow the definitions of the momentum 
scales, accepted in [18]: the hard scale, corresponding to momentum of order T, the soft scale ~ gT, and 
the ultrasoft scale ~ g 2 T. 
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equation for soft (ultrasoft) collective modes (in the case of Bose excitations - transverse 
and longitudinal modes, the latter are called by plasmons). Most efforts were directed 
at the derivation of the first type of equations. However for considerably excited states, 
when a characteristic time of relaxation of the hard particle distributions f s ,s = q,q,g 
is commensurable with a characteristic time of relaxation of the soft oscillations or even 
significantly exceeds it, along with kinetic equations for hard particles it is necessary to 
use the kinetic equation for the soft modes. 

The calculation of the collision term for the quark-gluon plasma was apparently first 
made in [11]. Within the framework of the concepts developed in the theory of electron-ion 
plasma, the scattering probability of (anti) quarks among themselves through longitudinal 
and transverse virtual oscillations which account for dynamical screening, is derived. We 
note that in this paper the Boltzmann equation for hard quarks and antiquarks was 
supplemented by an equation characterized the relaxation of field excitations and was 
presented as the second type of kinetic equation mentioned above. In paper [12] a similar 
Boltzmann equation for hard gluons was obtained. Here, the scattering probability was 
deduced within the framework of usual diagrammatic perturbation theory with inclusion 
of screening effects in the random-phase approximation (one-loop order). Although the 
relativistic Boltzmann equations constructed in these papers take into account such an 
important QGP property as screening, the range of their validity restricts their use to 
just colorless deviations from equilibrium distribution functions. 

In [4] within the framework of the (semi) classical representation of QGP [13] the 
Balescu-Lenard-type collision terms for small color and singlet deviations of the distri- 
butions from the initial colorless equilibrium, were derived by the Klimontovich method. 
However, in these papers, the organizing role of the various momentum scales was not 
recognized, resulting in some inconsistent and complicated transport equations for hard 
particles. 

Only in recent years it been possible to derive the Boltzmann equation for hard modes 
of hot non-Abelian plasma by a rigorous and self-consistent way using an expansion in the 
coupling constant, and clearly clarify the nature of the approximation involved, and thus 
fix its range of applicability. It was be shown that for longer wavelengths (A ~ l/g 2 T) 
of color excitations in the non-Abelian plasma, not only was consideration of interaction 
of hard particles with the soft degrees of freedom represented by mean fields essential, 
but also taking into account the collisions of hard particles among themselves. A similar 
Vlasov-Boltzmann equation reproduces exactly (at leading order in g) a large variety 
of the thermal results obtained by a more fundamental analysis of the diagrammatic 
perturbation theory [14] and provides in some cases a letter description of phenomena do 
not yield to a perturbative analysis. 

Here, one can extract three approaches to constructing an effective kinetic equation for 
hard particles with collision term. The first of them is connected with Bodeker's effective 
theory for the ultrarelativistic field modes [15]. Starting from the collisionless non-Abelian 
Vlasov equation, which is the result of integrating out the scale T [9], Bodeker has shown 
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how one can integrate out the scale gT in an expansion in the gauge coupling g. At leading 
order in g, he has obtained the linearized Vlasov-Boltzmann equation for the hard field 
modes, which besides a collision term also contains a Gaussian noise. Subsequently, this 
equation was also proposed by Arnold et al [16] who derive the relevant collision term on 
phenomenological grounds - by analyzing the scattering processes between hard particles 
in the plasma. The kinetic equation derived in [15, 16] has a non-trivial matrix structure, 
since, the distribution function that describes color fluctuations is not diagonal in color 
space. 

Afterwards, an alternative derivation of the collision term of Balescu-Lenard-type was 
proposed by Litim and Manuel, and by Valle Basagoiti [17]. The former authors used 
a classical transport theory, whereas the latter used the set of "microscopic" dynamical 
equations coming from the HTL effective action describing the evolution of the collision- 
less plasma. In both cases the collision terms were derived by averaging the statistical 
fluctuations in the plasma on the basis of the method developed by Klimontovich [3]. 

Blaizot and Iancu [18] suggested a detailed derivation of the Vlasov-Boltzmann equa- 
tion, starting from the Kadanoff-Baym equations. The derivation is based on the method 
of gauge covariant gradient expansion, which was first proposed by them for the collective 
dynamics at the scale gT [9]. By using the given equation they obtain [19] the effec- 
tive amplitudes for the ultrasoft color fields, which generalize the HTL's by including the 
effects of the collisions (see also Guerin [20]). 

The paper of Bezzerides and DuBois [21] devoted to the non-thermal QED plasma, is 
one of the first papers, in which the relativistic kinetic equation for the soft correlation 
function was considered. Over 60-70 years in connection with application to thermonu- 
clear fusion in the theory of the nonrelativistic electron-ion plasma, a powerful pertur- 
bative method (the so-called weak turbulent approximation) was developed [22-24] for 
research into various nonlinear plasma processes of the following types: on-shell scatter- 
ing of soft modes on hard particles (the nonlinear Landau damping), three- and four- wave 
decays, etc. In spite of the fact that this method is not able to describe the phenomena 
connected with strong turbulence; nevertheless, it enables one within the framework of a 
unified scheme to encompass a wide class of plasma phenomena. The paper [21] is the 
attempt at the extending of the above-mentioned weak turbulent theory to the electron- 
positron-photon plasma governed by the quantum electrodynamics. As a basic tool for 
such an extension, the CTP-formalism was used. However, since the main effort here 
was directed into investigating the collision integrals for hard electron and positrons, the 
authors have restricted their derivation to the plasmon kinetic equation, only taking into 
account pair production and the (linear) Landau damping process. Properly high-order 
processes, which we are interested in and are responsible for the nonlinear interaction 
mechanisms of the plasma waves, were not considered at all (see also the discussion of 
one-loop computations below). 

A similar kinetic equation for soft Bose-modes in a QGP was apparently first made by 
Heinz et al [25]. In the context of the imaginary time formalism in the one- loop approx- 
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imation the imaginary part of the complete color linear response function was deduced 
and it is shown that it can be expressed in the form of the Boltzmann-Nordheim collision 
term. On the basis of such a derived rate of decay (T^) and the rate for regeneration 
of the perturbations (Tj), the kinetic equation defining the evolution of a phase space 
distribution iV(x, t; k, uj) of soft electric perturbations of the momentum k = (u, k) in the 
form proposed by Weldon [26], was written as 

^ = -AT d + (i + iV)rv (l.i) 

at 

As in a previous case, the higher orders, when and Tj itself can functionally (in the 
general case, nonlinearly) depend on N here, were not considered. However the derivation 
of such dependence becomes important if we take into account that all computations in 
[25] ([21]) were performed with the "rigid" one-loop approximation, with bare propagators 
of massless gluons and quarks (electrons). However, as is known, quarks (electrons) and 
gluons inside the loop are not massless, they acquire the effective temperature-induced 
masses. The consequence of this fact is a kinematic prohibition of a decay of the soft 
perturbations into physical states. Furthermore by virtue of the fact that the phase 
velocities of both transverse and longitudinal eigenmodes of the plasma exceed velocity 
of light, the linear Landau damping is also absent. By virtue of the above-mentioned, the 
rates of decay and regeneration are just zero in this approximation. 

This paper is devoted to further study of the kinetic equation for soft modes of the 
non-Abelian plasma. The theoretical framework of this paper is derived from synthesis of 
two formal developments. The first one is the development of a nonlinear theory of plasma 
wave interactions in ordinary plasma - more exactly the weak turbulent approximation 
[22-24]. The second is the development of effective theory of hot QCD originally proposed 
by Braaten and Pisarski [27], Frenkel and Taylor [28], Jackiw and Nair [29], Blaizot 
and Iancu [9] and then developed in the papers [15-20]. In our previous papers [10] 
without resorting to a complicated diagrammatic technique within the framework of the 
semiclassical representations, the following term in the expansions of Td and Tj, linear 
on a phase space distribution of soft perturbations, was derived. However, as was shown 
in [10] (see also section 2), this approximation is not sufficient for a complete definition 
of the relaxation process of the soft modes in QGP. Here, we consider the next terms in 
the expansions of and Tj, and show that the corresponding nonlinear equation ( |1 . 1| ) 
is of purely Boltzmann type, i.e. the collision term on the r.h.s. of this equation has a 
standard Boltzmann structure, with a gain term and a loss term. 

The outline of the paper is as follows. In section 2 the preliminary comments, with re- 
gard to derivation of the Boltzmann equation, describing the plasmon-plasmon scattering 
are explained. In section 3 the essential features of the scheme, which we used previously 
in [10] to derive the kinetic equation with allowance for the nonlinear Landau damping, 
are summarized. In section 4 we discuss the consistency with gauge symmetry of the 
approximation scheme used. Section 5 is devoted to the determination of the interacting 
fields in the form of the expansion in free fields with the necessary accuracy for further 
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research. In section 6 we select all terms in the expansion of the color random current 
responsible for the four-plasmon decay and derive the intermediate kinetic equation which 
then in section 7 will be rewritten in the terms of HTL-amplitudes. Section 8 is devoted 
to deriving the probability of plasmon-plasmon scattering, which is the main result of this 
work. In the next section on the basis of the explicit form of the obtained collision inte- 
gral, the expression for lifetimes of colorless plasmons is defined and an estimate for the 
leading order in the coupling at the soft momentum scale is deduced. Finally in section 10 
we present our conclusions and future avenues of study. 



2 Preliminary comments 



We denote the localized number density of the plasmons by N l (k,x) = iV"£, and the 
distribution function of hard thermal gluons by f(p,x) = f p . In this paper we consider 
processes with longitudinal oscillations only, propagating in a purely gluonic plasma, with 
no quarks. Besides, we suppose that there is no external color current and/or mean color 
field in the system, and the system is in the global equilibrium state, i.e. 



tab 



ab . 



,E P /T _ I 



(2.1) 



Here. 



E p = |p| for a massless hard gluon, the coefficient 2 takes into account that the 



hard gluon has two helicity states and a,b — 1, . . . , N 2 — 1 for the SU(N C ) gauge group. 
The triviality of the color structure of the plasmon number density is a consequence of 
these restrictions (see section 4) 
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The dispersion relation for plasmons u = u (k) 



is defined from 



where 



(w, k) = 1 + 



3u. 
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F(x) = 



x 



In 



1 + ar 



x 



- m6(l - |x|) 



(2.2) 



(2.3) 



is longitudinal color permeability and Up t = g 2 N c T 2 /9 is a plasma frequency. 
We expect the time-space evolution of iV£ to be described by 
f)N l f)N l 



(2.4) 



where = du^/dh is the group velocity of the longitudinal oscillations^. For a general- 
ized decay rate and inverse decay rate Tj it is shown in an explicit form that in general 

2 Notice, that in the general case, when the distribution function of hard gluons is a slowly varying 
function i n ti me and space, equation ( |2.2| ) is replaced by Ree z (o;, k; t, x) = 0. In this case the l.h.s. of 



equation (2.4) should be supplemented by (duj 1 (k, x)/dx) (dN l (k, x)/dk) . 
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case they are functionals dependent on the plasmon number density. Although the ap- 
proach we shall use in the subsequent discussion is correct only within the framework of 
semiclassical approximation, it is convenient to interpret the terms entering into and 
r i; using a quantum language. 



The Eq. (|2.4|) in general, describes two principal processes of the nonlinear wave- 



interaction. The first of them represents the process of the stimulated emission and 
absorption of the collective wave quanta by hard particles of plasma. In this case the 
more general expression for the decay rate and the regeneration rate can be 
written in the following forms, respectively: 

rfM = E / tSs / dT - ™(pl k > k i> • • • > k - K, ■ ■ ■ , k'J < - - - Al (2.5) 



(27 

x(l + JV£,)... (l + #ij/p[l + /p 



and 



rPlNl) = W / ar(% w( P >\k[, . . . , k' m , k, k 1; . . . , k n ) jv£, ...Nl, m (2.6) 

n,m J \ Z7l J J 



x(l + Ar<J...(l + <J/ p ,[l + / ] 



PJ 



Here, the phase-space integration is 

J nrn ~ J (2vr) 3 ' ' ' (2tt) 3 (2vr) 3 ' ' ' (2tt) 3 { } 

x (2tt) 5(£p + ^ + < + ... + < - # P ' - 4; " • • • " 4j , 
with the delta function expressing the energy conservation of the processes of stimulated 
emission and absorption of the plasmons. The function w(p|k, k 1; . . . , k n ; k' 1; . . . , k^) is 
the probability of absorption of n + 1 plasmons with the frequencies u> l k , uj{. v . . . , uj l kn and 
the wavevectors k, ki, . . . , k n by a hard gluon carrying of momentum p with consequent 
radiation of m plasmons with frequencies u^, , . . . , and the wavevectors k' 1; . . . , k' m . 
The function u^p'm, . . . , k' m ; k, ki, . . . , k„) is the probability of inverse process - the 
absorption of m plasmons by a hard gluon with the momentum p' = p + k + ki + ...+ 
k n — k' x — . . . — k.' m with consequent radiation of n + 1 plasmons. Diagrammatically this 
corresponds to a Feynman graph with two hard external lines and an arbitrary number of 
n + m + 1 soft external lines. By using the fact that |p| ^> |k|, |ki|, . . . , |k n |, |m, . . . , \k' m \, 
the energy conservation law can be represented in the form of following "generalized" 
resonance condition: 

4 + ^ i +...+4 n -^,-...-^-v(k + k 1 + ... + k n -k' 1 -...-k^) = 0, (2.8) 

where v = p/|p| . Furthermore, one can approximate the distribution function of hard 
gluons on the r.h.s. of equations ( |2.5| ) and (|2.6|), 

df 

/ p ,~/ p + (k + k 1 + ... + k n -k' 1 -...-k'J^, 
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and set 1 + f p ~ 1 + f p > ~ 1 by virtue of f p , f p > <C 1. 

The r.h.s. of (|2.5|) and (|2.6|) can be formally considered as an expansions of 
and in the functional series in powers of the plasmon number density. The actual 
dimensionless parameter of expansion here, is (for classical statistic) the ratio of the energy 
of longitudinal plasma excitations to the averaged thermal energy per particle, i.e. 



where n is the mean density. In conditions, when the excitations energy is a small quantity 
compared with the thermal energy of hard particles, we have 

£<1. (2.9) 

The last inequality means that the fields of longitudinal oscillations are sufficiently small 
and they cannot essentially change such "crude" equilibrium parameters of a plasma as 
particles density, temperature and thermal energy (this, in particular, justifies the choice 
of the distribution function of thermal gluons in the form of ( |2.1| )). 

On the other hand, however, we shall consider the energy of the plasma oscillations to 
be sufficiently large, i.e. greatly exceeding the energy of thermal fluctuations of the the 
color field in the plasma. The consequence of the last requirement is the inequality 



e > S, (2.10) 



where 5 is the plasma parameter 



^.3 

5 = ^- < 1. 

r 6 D 



Here, r is the inter-particle distance (~ n 1//3 ), and is the Debye length 

T 



r D 



4nng 2 N c 



The condition (|2.10|) demonstrates the validity of ignoring the hard gluon collisions among 



themselves relative to their interactions with soft plasma modes. 



Inequalities fl2.9|) and (|2.10| ) correspond to the weak turbulent approximation, within 



the framework of which one can restrict the consideration to several first terms in a 

(s) (S) 

functional expansions of T A and V\ . We note however, that when the energy level 
of the plasma excitations becomes comparable with the thermal energy of particles, e.g. 
as the result of development of a strong instability (strong turbulence), the perturbation 
theory here, is no longer applicable and the problem of summation of all the relevant 
contributions thus appears. The last situation can be really take place in the processes 
proceed in QGP emerging from the heavy ion collisions at higher energies. In this work 
we do not consider this very complicated problem, assuming that the inequality ( |2.9| ) is 
always fulfilled. 
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Let us discuss in more detail the first terms in the expansions of r£ ) and rf } . For 



n = m = 0, the equation (|2.8f) results in the relation 

<4 - vk = 0, 

which is well-known as £/je Cherenkov resonance condition, which does not hold in the 
gluon plasma. Therefore Tffl = v\ = 0, when only 0(e°) terms on the r.h.s. of 
and (12.0) are retained. 



For n — 1, m — we have 

4 + <-v(k + k0 = 0, (2.11) 

and for n = 0, m = 1, respectively (here, we replace k[ by ki) 

4-< -v(k-ki) = 0. (2.12) 

The first resonance condition (|2.11|) describes the simultaneous radiation (or absorption) 
of two plasmons with frequencies u^, and wavevectors k, ki. By virtue of the fact 
that the phase velocity of the longitudinal oscillations exceeds the velocity of light, this 
process is kinematically forbidden. The first non-trivial terms in the expansion of the 
functionals rf and rf ^ are defined by the second resonance condition (|2.12| ). It is 
associated with the absorption of the plasmon by a hard gluon with frequency u l k and 
wavevector k with its consequent radiation with frequency uj^ and wavevector ki (and 
vice versa). Schematically this process can be represented as follows: 

g* + g^g*+g, (2.13) 

where g* are the plasmon collective excitations and g are excitations with characteristic 
momenta of order T. In the theory of the ordinary plasma [22-24] this process is known 
as the nonlinear Landau damping. In the case of QGP it was studied in detail in [10]. We 
have shown, that the nonlinear Landau damping rate 

V(k) = (rf[Nl] - rf \Nl])\ n=0 , m=1 (2.14) 

defines two processes: the effective pumping of energy across the spectrum towards small 
wavenumbers with the conservation of excitation energy and properly nonlinear dissipa- 
tion (damping) of the longitudinal plasma waves by hard particles, where the first process 
is crucial. The consequence of this fact is the inequality: 7*(0) < 0, i.e. k = - mode 
is increased. The main conclusion, which we drew in [10] is that the only process of the 
nonlinear Landau damping does not lead to the total relaxation of soft excitations in the 
gomogeneous isotropic plasma. At the scale of a small |k| (|k| <C gT) it is necessary to 
consider the processes of higher-order in e, than ( 2.13Q , which lead to the suppression of 
increase of the k = 0-mode. 

The following terms in the expansion of the decay rate ( |2.5|) and the inverse decay 
rate ( ^.6| ), corresponding to n,m = 1,2, are defined by 

4 ± ui ± uL - v(k ± kx ± k 2 ) = 0, (2.15) 
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± 4 =F < - v(k ± ki T k 2 ) = 0. 



Physically this corresponds to simultaneous absorption (radiation) of three plasmons by 
a thermal gluon, or simultaneous absorption (radiation) of two plasmons with consequent 
radiation (absorption) of one plasmon. At the long-wavelength range these processes are 
kinematically forbidden and therefore in this approximation and r[ S ' vanish. 

However there are contributions different from zero in the expansions of the general- 



ized decay rate and the inverse I\ in the equation (2.4) which are of the same order 



as the processes Q2.15Q , i.e. of order 0(e 2 ). These contributions are concerned with the 



second type of the nonlinear processes defining the time-space evolution of iV£ and going 
without exchange of energy between hard thermal gluons and plasmons. They represent 
the processes of decays, fusions of plasmons and their scattering off each other. Dia- 
grammatically, this corresponds to Feynman graph, where all external lines are soft. The 
relevant decay rate Y\ and regenerating rate T d can be formally represented in the 
form 

rf } [iV£] = £ / dTffl w(k, k 1; . . . , k n ; k' x , . . . XX . . . Nl n (l + Nl,) . . . (1 + N&, 



n,m ' 



(2.16) 
and 

r?V£] = E / dT Z } w(K, . . . , k' m , k, k x , . . . , k„) Nl, . . . Nlji + Nlj . . . (i + Nlj. 

n,m 

(2.17) 

Here, the phase-space measure is 

(2.18) 

with the delta functions expressing the energy and the momentum conservation of the 
decay processes, the fusion and the plasmon scattering among themselves. The decay rate 
( p,16| ) and the regenerating rate (|2.17| ) are different from zero if the following " resonance" 
conditions are obeyed 



uL + uL+... + ui-uL-...- uL = 0, 



k + ki + . . . + k n - k' x - . . . - k^ = 0. 

The first contribution is different from zero in and , which would arise in the 
case of n = m = 1, 2 as defined by the system of equations 

4±<T<=0, (2.19) 
k ± ki =f k 2 = 0. 

These conservation laws describe a decay of one plasmon into two plasmons and the 
reverse process of fusion of two plasmons into one plasmon, 

g*^gl + g* 2 - (2.20) 
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Three-plasmon decay ( |2.20| ) is as important as the process of nonlinear Landau damping 
( [2.13] ). However, the specific peculiarity of a dispersion law of the longitudinal oscillations 
in hot non-Abelian plasma is that resonance equations ( |2.19p have no solutions no matter 



what the values of wavevectors k, ki and k 2 may be. The processes ( |2.20|) are kinematically 
forbidden by the conservation laws, and it makes a contribution only in the second order 
over parameter e of perturbation theory, for n, m — 1, 2, 3. 

In general case, with four plasmons, two different processes occur 

g*#gl + g 2 * + g^, (2.21) 

g* + g^ g ; + g*. (2.22) 

The first of them corresponds to the process of the decay of one plasmon g* into three 
plasmons g*, g 2 , g3, and the reverse process of fusion of three plasmons into one plasmon 
g*. The second process presents the plasmon scattering by plasmon. The last one is 
considered as the decay process and is interpretated as the process of the decay (fusion) 
of two plasmons g* and g^ into two plasmons g 2 and gg. For four-plasmon decay process 
( p.21| ), the following resonance conditions are obeyed: 

4-<-<-< = 0, (2.23) 

k - ki - k 2 - k 3 = 0, 
and for the scattering process (|2.22|) we have, respectively, 

4 + <-4 2 -4 3 = 0, (2.24) 
k + k x - k 2 - k 3 = 0. 



It is not difficult to show that the conservation laws fl2.23|) and (|2.24|) kinematically forbid 



the processes of the direct decay of the plasmon into three (and vice versa) Q2.21| ) and 
admit only the processes of ( |2.22|) type. As shown in section 9, at the soft scale the latter 
process is suppressed by a power of g relative to the process of nonlinear Landau damping 
( p. 13 ). However at the ultrasoft scale, one would expect that the process of the elastic 



scattering of the plasmon by a plasmon may be as larger as the process ( |2.13| ) and thus 



plays an important role in the kinetics of the plasmons at larger wavelength (see the end 
of section 9). 



Putting the expressions ( |2.16|) and (|2.17|) into equation (|2.4|) , where only 0(e 2 ) relevant 



terms in the expansions are retained, we result in the Boltzmann equation, describing 
four-plasmon decays of ( [2.22j ) type 

dNL , dNl . . 

-§± + VL 1 £= (2-25) 
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xw(k, ki; k 2 , k 3 ){«(l + Nl)(l + JVjJ - iVX^l + + O) 

(27r)45(4 + ^ " 42 " ^ )5(k + kl " k2 " k3) 

x W (k, k i; k 2 , kB){«^i, + - NlNlNl 2 - NlNlNij. 

In writing this equation, we have used the fact that the probabilities of direct and reverse 
processes are equal and besides in the last line in the semiclassical regime we consider the 
soft modes to be strongly populated, i.e. (1 + iV£)(l + A^J « N^N^ +N^ + etc.. A 
similar Boltzmann equation for plasmons was studied intensively for the ordinary plasma 
[30, 31] (in [24, 30] the explicit expression for the function w(k, k^ k 2 , k 3 ) can be found). 
The main purpose of this work is the derivation in the explicit form of the probability of 
plasmon-plasmon scattering for hot non-Abelian plasma. The function of w(k, k x ; k 2 , k 3 ) 
must satisfy the symmetry relations over permutation of arguments 

iw(k, ki; k 2 , k 3 ) = w(k 2 , k 3 ; k k x ) = w(k, ki; k 3 , k 2 ) = w(ki, k; k 2 , k 3 ), (2.26) 

which are the consequence of the indistinguishable of the plasmons (recall that here, we 
discuss only colorless plasmons, i.e. those for wich the occupation number does not carry 
adjoint color indices). 

In conclusion of this section we note some general properties of the processes of four- 
plasmon decays. Multiplying the r.h.s. of Eq. (2.25[) in turn by uj\, and k, integrating 



with respect to the wavevector k, and taking into account ( 2.26Q , it is easily checked that 



/dk r dk 

— J k N^ const, K .y_k< = const. 

These relations are evident a consequence of the conservation laws of energy and momen- 
tum in the plasmon-plasmon scattering. 

On the other hand, the total plasmon numbers in such decays should also be conserved, 
because the only process in which two plasmons decay into two others is permitted. Really, 
by integrating Q2.25D over all k-space and taking into account the relations fl2.26p , it is 
easy to verify that 

i.e. four-plasmon decay not changes the total number of plasmons. 



3 The random-phase approximation 

In this section let us briefly recall the main methodological notion used in study of the 
processes of nonlinear wave-interaction in a non-Abelian plasma within semiclassical ap- 
proximation in [10]. 
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We use the metric g^ v = diag(l, —1,-1,-1) and choose units such that c = ks — 1. 
The gauge field potentials are N c x A^-matrices in a color space defined by = A a ^t a 
with N^ — l Hermitian generators of the SU (JV C ) group in the fundamental representation. 
The field strength tensor = F^t a with 

F; v = d»K - d v Al + gj ahc A\Al (3.1) 
obeys the Yang-Mills (YM) equation in a covariant gauge 

dpF^X) - ig[A^(X), F^{X)\ - C^d'A^X) = -f(X), (3.2) 
where £ is a gauge parameter. j u is a color current 

f = gt a J d 4 pp»Tr(T a f), (3.3) 

where T a are Hermitian generators of SU (N c ) group in the adjoint representation ((T a ) bc = 
-if abc ,Tr(T a T b ) = N c 5 ab ). The distribution function of gluons / satisfies the dynamical 
equation which in the semiclassical limit (when polarization effects are neglected [32]) is 
[13] 

P^m/+^{^>.^} = 0» ( 3 - 4 ) 

where is a covariant derivative acting as 

V ll = d ll -ig[A ll (X), •], 

with [ , ] and {, } denoting the commutator and anticommutator, respectively, and A^, 
J> are defined as A» = A*T a , = F^T a . 

The distribution function / can be decomposed into two parts: regular and random, 
where the latter is generated by spontaneous fluctuations in the plasma 

f = f R + f, (3.5) 

so that 

(/> = f R , if) = o. 

Here, angular brackets (•) indicate a statistical ensemble of averaging. The initial values of 
parameters which characterize the collective degree of plasma freedom is a such statistical 
ensemble. For almost linear collective motion to be considered below this may be initial 
values of oscillation phases. 

We also use the definition 

A, = + Al , {AD = 0. (3.6) 

The regular (background) part of the field A^ will be considered to be equal to zero. 
The condition for which the last assumption holds, will be considered closely in the next 
section. 
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By averaging equation ( |3.4j ) over statistical ensemble, we obtain the kinetic equation 
for the regular part of the distribution function of hard gluons / 

1 8f T 1 8f R 

P^f R = igiT([Al, f]) - 2^({(^)x, ^-}> - ^{((J^W, 

(3-7) 

Here, the indices L and iVX denote the linear and nonlinear parts with respect to field 
A a of the strength tensor ( |3.1| ). The correlation functions on the r.h.s. of this equation 
are collision terms due to particle-wave interactions and describe the backreaction of the 
background state from the plasma waves. 

We assume that the typical time the nonlinear relaxation for the oscillations is a small 
quantity relative to the time scale over which the distribution of hard transverse gluons f R 
vary substantially. Therefore we neglect by change of the regular part of the distribution 
function with space and time, assuming that this function is specified and describes the 
global equilibrium in the gluon plasma 

/W = 2 |M V) _±_ (3 . 8) 

where is the 4-velocity of the plasma. (Here, for convenience, we somewhat overdeter- 
mine the equilibrium distribution function of thermal gluons (|2.1| ).) 



We use the expansion in powers of the oscillations amplitude of the random function 
f T to investigate non-equilibrium processes in QGP, such that the excitation energy of 
waves is small quantity in relation to the total energy of the particles 



f=Y.f {n \ (3-9) 



pT _ f T(n) 
n=l 

where f T ^ collects the contributions of the ra-th power in A T . The expansion of a color 
current, corresponding to (|3.9|) has the form 



U = i? + il , <j M > = $ , jI = E ij (n) > (3-io) 

n=l 

where by the definition ( |3.3|) , we have 

ij (n) = 9t a J d 4 pp,Tr(T a f T ^). (3.11) 
The regular part of a current vanishes for the global equilibrium gluon plasma. 



Substituting the expansion ( |3.9| ) into (|3.4|) , and collecting the terms of the same order 
in A T , we derive the system of equations 

^-^{W^L (3-12) 
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P"d,f^ = igjf([Al, fM] - ([A*, f T ^]}) (3.13) 
P^f^ = igpi*([Al, / T(2) ] - ([Al / T(2) ]» - \g P ^m^ ^} (3.14) 



We rewrite the Yang-Mills equation ( |3.2| ), connecting a gauge field with a color current, 
in the following form 

d^F T ^) L - C^d^Al + j 2 ™* = (3.15) 

= -J T nl + i9d»([A T », A Tv ] - ([A T », A Tv })) 

+W([AI, (F T n L ] - ([Al (F T n L ]))+9 2 ([Al, [A T »,A**\] - (L4j, [A^,A T »]})). 

Here, on the l.h.s. we collect all linear terms with respect to A^ and we denote: jJ^ L = 
jT(2)v _j_ jT(3)u _)_____ w jii shown that at leading order in g only, the first two terms 
in jJf U L need to be kept. 

It is not difficult to obtain the explicit form of the terms in the color current expansion 
( p.lOj ) from the system of equations ( |3.12j )- (|3.14 ) and the relation ( 3.11 ). In the momen- 
tum representation and to the leading order in coupling constant up to the third-order 
terms, we have [10] 

f^(k) = U^(k)A T u (k), (3.16) 

where 

n~(*) = f 6^M, e +0 (3.17) 
J pk + ip e 

is the high temperature polarization tensor; 

fW*(k) = f abc J &&£{Jb{k x )JTM (3- 18 ) 
where 

s£2T = -V". i ** (3-19) 

' u 2 J pk + ip Q e pk 2 + ipoe 

f&X*( k ) = r *f f f*eJ ECJ2Xi(^(*»)^(*i)^(*2) - 4(^)(^(A;i)^(A; 2 )) 

- (^(A; 3 )^(A; 1 )^(A; 2 )))5(A; - A* - A; 2 - fe 3 ) dhdk 2 dk 3 , (3.20) 

and 

= 4jV / rf 4 P^pY 1 (W°) (3 21) 

For simplicity, hereafter we drop the superscript T of a gauge field. 
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Furthermore, we rewrite the equation ( |3.15| ) in the momentum representation. By 
inserting the linear part of the random current ( |3.16| ) and nonlinear corrections (|3. 18 ) 
and ( |3.2(J| ) into Eq. ( |3.15| ), one finds 

[k 2 g^ - (1 + C 1 )^^ - W v (k)}A b v (k) 

= f cd J S£t M {Al^) Ad x^2) - {AUkMli^mk -h- k 2 )dhdk 2 (3.22) 

+f bcf f fde J E^,J^(A*)A*(A^ 

x5(k — k\ — k 2 — k%) dk\dk 2 dk^. 

Hprp = Q( / )^ 1/A i c( /7 W yfJ.f^ y(I)iiv\a y (II)nv\a rp, f 11T , rt ; nT ,„ q(II) 

neie, o k k ^ k2 _ o k ki ^ k2 -|- o k kl k2 , ^k MMM — ^kMMM * ^kMMM' e luncuons o 

and X^ 7 ) are defined by expressions (|3.19| ) and Q3.21|) , respectively, and 

S%!£ = -ig(Wg* + - fc 2 V A ), ^ttM = sVV* (3.23) 

These tensor structures are caused by self-action of a gauge field. They are defined by 
nonlinear terms on the r.h.s. of equation ( 3.15 ), which are not associated with a color 
current in QGP. In section 5 equation Q3.22j ) will be formally solved by iteration. 

In [10] we introduce the correlation function of random oscillations 

l£(k',k) = (Al a (k')Al(k)). (3.24) 

In order not to over burden equations by the symbol " * " we use a dagger " f " to denote 
complex conjugation. In thermal equilibrium, when the correlation function ( |3.24| ) in 
the coordinate representation depends only on the relative coordinates and time AX = 
X' — X, we have 

I^(k',k) = I^(k')5(k'-k). (3.25) 

Off-equilibrium perturbations which are slowly varying in space and time lead to a 
delta function broadering, and 1°^ depends on both arguments k and k 1 . 

Let us introduce I^,(k',k) = I^l(k, Ak), Ak = k' — k with | Ak/k |<C 1 and insert 
the spectral intensity function in the Wigner form 

l£(k,x) = J I^ b u {k, Ak)e~ iAkx dAk, 

depending slowly on x. 

Now we multiple equation ( p.22| ) by A^(k'), subtract to it the complex-conjugated 
equation (with the replacement k *->■ k' , a b) and average the equation using the 
formula ( |3.24j ). Furthermore we expand the polarization tensor into Hermitian and anti- 
Hermitian parts 

n w (fc) = ii Hva {k) + ii AycT (k) , n Hva (k) = Yi ]Hav (k) , ii Ava (k) = -itf Aav (k). 
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The term with U A corresponds to linear Landau damping. As was shown by Heinz 
and Siemens [33], that linear Landau damping for waves in QGP is absent and hence, 
this term vanishes. We expand the remaining terms on the l.h.s. in powers of Ak and 
keep only linear ones. This corresponds to agradient expansion procedure, usually used 
in the derivation of kinetic equations. Multiplying the difference equation by e~ lAkx and 
integrating over Ak with regard to 

J Ak\ Zg(fc, Ak) e-^dAk = z 91 ^ X) , 

we finally obtain the equation, which is a starting point for our further research 

r) dl ab 

[k 2 g^- (l + r 1 )*^ -H H ^(k)} — ^ = (3.26) 



dk\ x dx x 

= -%J dk'{f^ d S& lM (Al a (k')Al(k^ 

xS(k — k\ — ki — k^dk\dkidk?, — (a <-> b, k <-> k', compl. conj.)}. 



4 Consistency with gauge symmetry 

In this section we shall discuss the consistency of the approximation scheme which we use 
with the requirement of the non-Abelian gauge symmetry. 

The initial dynamical equation ( |3.4| ) and the Yang-Mills equation (|3.2| ) (without the 
gauge-fixing condition) transform covariantly under local transformations 

A^X) = h(X)(A fl (X) + -djrfiX), h(X) = exp {i9 a {X)t a ) 

with the parameter 6 a (X). We also have transformation of gluon distribution function 
[13] 

f(p,X) = H(X)f(p,X)H\X), 
where H ab (X) = Sp[t a h(X)t b tf(X)}. 

As is known (see, e.g. [17]), after the splitting ( |3.5|) , ( |3.6|) the resulting equations left 
two symmetries: the background gauge symmetry, 

A$(X) = h(X)(A*(X) + ~d„)h){X), AJ(X) = h{X)Al{X)h\X), (4.1) 

and the fluctuation gauge symmetry, 

A«(X) = 0, AJ(X) = h(X)(A*(X) + Al(X) + - g d,)h\X). (4.2) 
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The condition which we impose on a regular part of the gauge field in the preceding 
section and the requirement that the statistical average of the fluctuation vanishes (A^) = 
0, break down both types of symmetry ( [4.1| ) and Q4.2p. Thus in the case of a gauge 
transformation ( |4.1|) we obtain A^ ^ 0, and in the case of ( |4.2| ) we arrive at non-invariance 
of the constraint {AV) = 0. Moreover, the introduced correlation function (|3.24|) also has 



an explicitly gauge non-covariant character. This leads to the fact that calculations in the 
previous section are gauge non-covariant, and therefore the value of these manipulations 
is doubtful. 

Nevertheless, there is a special case, when the preceding (and following) conclusions 
are justified. This is a case of colorless fluctuation, where I®l(k,x) = S^I^k^x). We 
can obtain a gauge-invariant equation for I^kjx) only in this restriction, in spite of the 
fact that the intermediate calculations spoil non-Abelian gauge symmetry of the initial 
equations (|3.2|) 



In principle, we shall be able to maintain an explicit background gauge symmetry 
Q4.1D at each step of our calculations, as has been done, for example, by Blaizot and Iancu 
[18] for derivation of the Boltzmann equation describing the relaxation of ultrasoft color 
excitations. First of all we assume that A^ ^ 0. Then as the gauge-fixing condition for 
the random field A T , we choose the background field gauge 

V*{X)A T »{X) = 0, V*{X) =d„- igA*(X), (4.3) 

which is manifestly covariant with respect to gauge transformations of the background 
gauge field A^(X). Lastly we define a gauge-covariant Wigner function as in [13, 18] 

ft(k,x) = J f£(s,x)e iks ds, s = X 1 -X 2 ,x=^(X 1 +X 2 ), 



where 



If v {s, x) = U aa ' (x, x + J) I a ' b \x + J, x - -) U b '\x - J, x) 



instead of the usual Wigner function I°£(k,x), whose 'poor' transformation properties 
follow from the initial definition lf v {X x ,X 2 ) = (Al a (Xi)A™(X 2 )). The function U(x,y) 
is the non-Abelian parallel transporter 



U(x,y)=Pexp{-ig J dz»A«(z)}. 



The path 7 is the straight line joining x and y. 

The derivation of the kinetic equation for plasmons in this approach becomes quite 
cumbersome and non-trivial. For example, on the l.h.s. of the equations for a random 
part of the distribution ( 3.12|) -( 3.14j) , the covariant derivative T>^ will be used instead of 



the ordinary one d^. Besides, we cannot assume that the regular part of the distribution 
function is specified and equal to the Bose-Einstein distribution ( |3.8| ). It is necessary to 
also take into account their change using the kinetic equation ( p.7|) with the collision terms 
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on the r.h.s. of ( |3.7| ), which describes the backreaction of the background distributions 



from the soft fluctuations. The correlators on the r.h.s. of equation (3/7) can be expressed 



in terms of the function J~J and the distribution of hard transverse gluons f R (p, X) only. 

However, if we restrict our consideration to the study of colourless excitations and 
replace the distribution function of hard gluons by its equilibrium value (|3.8j), then this 
leads to an effective vanishing of terms with mean field A R . This follows, for example, 
from the analysis of the derivation of the Boltzmann equation by Blaizot and Iancu 
[18]. Therefore, the simplest way to derive the kinetic equations for soft colorless QGP 
excitations is to assume A R = and to use the prime gauge non-covariant correlator ( |3.24D . 
In this case the background field gauge (fO| ) is reduced to a covariant one. The resulting 
Boltzmann equation for colorless plasmons will be gauge invariant if all contributions to 
the probability of plasmon-plasmon scattering at the leading order in g are taken into 
account (see also discussion in conclusion). 



5 The interacting fields as the functions of free fields 

Let us define approximate solution of equation ( |3.22| ) accurate up to third-order in the 
oscillations amplitude of the free field. For this purpose, it is convenient to write this 
equation in a more compact form. 

We introduce the following notation 
jT { 2 W{k) = fabc J S^iAlikJAlih) - (A b u (k 1 )A c x (k 2 )))5{k -k x - k 2 )dhdk 2 , (5.1) 

jT(3W (k) s fab f ffd eJ v^ M (Al(kMx(h)Al(k 2 ) - AKhXAHh^M) 

- (AKk^AKkx^MWik -h-k 2 - h) dk l( lk 2 dk 3 . (5.2) 

The expressions (|5.1|) and (|5.2j) present the nonlinear color currents, including the self- 
action effects of gauge fields, in contrast to (|3.18|) and ( |3.20|) . 

Using *V^ v {k) we denote the medium modified (retarded) gluon propagator, which in 
a covariant gauge has a form 

*V^{k) = -P^k) *A t {k) - Q^(k) *A l (k) + £D^(k)A°(k), (5.3) 

where *A t ' l (k) = 1/(A; 2 -IT''(A;)), U^k) = \W v (k)P^(k), U l (k) = W \k)Q „ v (k); A°(k) = 
1/k 2 . The Lorentz matrices in (|5.3| ) are members of the basis 

Pfj,v(k) = g^u — D^ u (k) — Q^ u (k), Q^ u (k) = ^ _ 2 /> \ > C^ u (k) = — - — , 

u \ K ) J —2k 2 u 2 (k) 

D^u = k^/k 2 , u^(k) = Pu^ - k^ku). (5.4) 
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Let us assume that we are in the rest frame of a heat bath, so that = (1, 0, 0, 0). 

Using the above introduced functions, the equation (|3.22|) can be rewritten in the form 

*V~ llxv {k)A a v (k) = -J T{2)a ^(A, A) - J T(3) ^(A, A, A). (5.5) 

The nonlinear integral equation (|5.5| ) is solved by the approximation scheme method - 
theweak free field expansion (small perturbations). Discarding the nonlinear terms in A 
on the r.h.s. of equation ( |5.5|) , we obtain in the first approximation 

*T>~ ltJ,v (k)Al(k) = 0. 

The solution of this equation, which we denote by A^ a (k) is the solution for free fields. 

Further keeping the term, quadratic in the field on the r.h.s. of Eq. ( |5.5| ), we derive 
the equation 

*V- ltxv {k)A a v {k) = -J T ^(A^,A^), 

where on the r.h.s. we substitute free fields instead of interacting ones. The general 
solution of the last equation can be given in the form 

A a ^k) = Af a (k) - *V^ u (k)J T{2)av (A^\A^). 

This approximate solution was used in research of nonlinear plasmon damping in QGP 
[10]. 

The following term in the expansion of interacting fields is defined from equation 
*V- l ^{k)A a v {k) = _j^(2)^(_*pjT(2) (A (o) } A^),A^)-J T ^(A^,-*VJ T ^ 2 \A^,A^)) 

-J T ^(A(°) ) # ) #). (5.6) 

Using the explicit expressions for the currents ( |5.1| ) and ( |5.2| ) , after cumbersome algebraic 
transformations, we obtain the form of the interacting field from the equation (|5.6|) with 
the accuracy required for our further calculations 

A*(k) = Af a (k) - *V^(k)J T{2)au (A {0 \ A {0) ) - *V^(k)J T{3)au (A^\A^\A {0) ). (5.7) 

Here, the third-order color current on the r.h.s. is defined by the expression 

^(3)^(0)^(0)^(0)) _ f a bfffde J S^ 2ifc3 (Af b (fc 3 )4 0)d (fcl)4° )e (^)- 

- A^\k 3 )(A^ d {k 1 )Af>{k 2 )))S(k -h-k 2 - h) dhdhdh, (5.8) 

where 

^kMMM — ^kMMM 2 7 ^ ^ D k,k 3 M+k2 D k,ki+kiM> Wfci+fe,fei,fe D fci+fc 2 ,fc2,fcJ> 

(5.9) 

and we take into account, that the third-order correlation function (A^°' A^- ' A^ ') vanishes 
by virtue of the fact that A^ represents the amplitude fully non-correlative gauge fields. 
The factor of |, in front of the second term on the r.h.s. of (|5.9|) arises from symmetrization 
with respect to permutation of the potentials A^> d {k\) and A^ e {k2) in the expression 
Q5.8| ). The current ( |5.8| ) may be interpreted as a certain third-order effective color current, 
in contrast to the initial 'bare' expression ( [5.2|) . 
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6 The correspondence principle 



For the determination of the probability of plasmon-plasmon scattering in a gluon plasma 
the method developed in the theory of the nonlinear processes in electron-ion plasma and 
known as the correspondence principle [24, 30], is usable. For the non-Abelian plasma this 
approach is especially effective in the temporal gauge, when we have closer correspondence 
with the electrodynamics of an ordinary plasma. The gist of this method is as follows. 

The change in the plasmon numbers, caused by spontaneous processes of four-plasmon 
decays only, is 

B + v ! MV P = rdk^dk^dk^ 

\ dt k <9x ) i J (2tt) 3 (2vr) 3 (2tt) 3 1 ' 

x (2tt) 4 5(^ + < - < - 4 3 )5(k + k! - k 2 - k 3 ) w(k, k i; k 2 , k 3 )<iV^iV^. 

This equation follows from equation (|2.25| ) in the limit of small intensity — *> 0. In this 
case the change of energy of the longitudinal excitations is 

dk dki dk 2 dk 3 




(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 



(6.2) 



x (2vr) 4 ( 5(4 + < - 4 2 - 4 3 )5(k -k,-k 2 - kg) uj l k w(k, k l5 k 2 , k^N^N^N^. 

On the other hand the value (d£/dt) s / represents the emitted radiant power of the longi- 
tudinal waves I 1 , which in turn is equal to the work done by the radiation field with the 
color current, creating it, in unit time 

l l = J dx <E'(x, t) J»(x, t)) = J (E^ Jl^)e-W-»* = (6.3) 

= 1 r da; do/ dk f 1 W taj ■ i(a/ _^ 

~2J 2tt 2tt (2tt) 3 k) ^( w ,k)/ k 2 V k ' w k ' w ' 7 

Here, E ai (x.,t) = —dA ai {yL,t)/dt is chromoelectric field in the temporal gauge. The sign 
on the r.h.s. of ( |6 . 3| ) corresponds to the choice of sign in front of the current in the Yang- 
Mills equation (|3.2|) . In conclusion of the last line (|6.3| ) we take into account that the 
Fourier- component of a field E£ = kE£ /|k| is associated with J£ w by the Yang-Mills 
equation 

k ' w ue l {u,k) |k| ' 

In order to define the probability of the four-plasmon decays, the correlation function on 
the r.h.s. of equation ( |6.3| ) has to contain terms of six-order in the free field A^°\ The 
required sixth-order correlator yields the color current J T ( 3 ) m ( |5.2| ) (more precisely, its 
expression in the temporal gauge). However here, it is also necessary to take into account 
the effects, that arise from iteration of the current J T ( 2 ) a * ( |5.1| ). Defining in this way all 
necessary contributions, making the correlation decoupling of the sixth-order correlators 
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in terms of pairs and next expressing next (A^A^) in terms of N l , we obtain an emit- 
ted radiant power I 1 . Comparing X 1 with (|6.2|) , one identifies the required probability 
w(k,ki; k 2 ,k 3 ). 

However this method encountered certain difficulties in deciding on the other gauges, 
e.g. the covariant gauge. Here, it is convenient for the definition of the probability of 
plasmon-plasmon scattering to start immediately from equation ( |3.26|) . Using the above- 
obtained expression ( |5. 7| ) for the potentials of interacting fields, by a simple search we 
extract all the sixth-order correlators responsible for four-plasmon decays of the type 
( p.22j ). The rule used to choose the relevant terms is defined by the simple fact that when 
we make the correlation decoupling of the sixth-order correlators in the terms of the pairs 
(in order to define the product N^N^N^) as a factor, the ^-functions arise in the form 

5U + ul - ui - ul)5(k + k! - k 2 - k 3 ). (6.4) 



Also the coefficient function preceding NLN^ N^ must satisfy properties ( |2.26| ). As 



will be shown below, these conditions are sufficient to calculate the plasmon-plasmon 
scattering probability in the covariant gauge (this rather cumbersome and physically not 
quite transparent approach is suited for other gauges). 

At first, we consider the contribution to the r.h.s. of the initial equation ( |3.26| ), asso- 
ciated with S-functions. Here, for convenience of further reference we write it separately 

x 5(k — hi — hi — k^dkidkidkz — (a <-> b, k *-> k', compl. conj.)}. (6-5) 

One can obtain the sixth-order terms in the free field by two waves. The first one is as 
follows. We substitute in turn the expression from the r.h.s. of ( |5.7|) , which contains only 
cubic terms in the potentials of free fields, instead of each potential of the interacting 
ones, i.e. 

A»(jfe) - -*V^(k)J T ^'(A^,A^,A^). 

We replace the remaining potentials by the rule v4£(fc 3 ) — > A^ c (k 3 ), etc. The second way 
is to substitute the quadratic term in A^ from the r.h.s. of ( |5.7| ) instead of any two 
potentials of the interacting fields, i.e. 

A*(k) -> -*D^(A;)J T(2) ^'(A (0) ,A (0) ). 

We replace the remaining potentials by free ones. It is necessary to look at all possible 
substitutions in both first and second ways. 

The number of terms appearing can be cut if we note, that it is need to keep only such 
terms in intermediate expressions which contain the propagators *T>^^(k) and *T> tlfl i(k'). 
These propagators give the terms, proportional to 

5(Ree\k)) = ( dRe d ^ ik) Y 1 , [5(u - ufo + % + 4)], (6-6) 
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i.e. the factor taken into account the existence of plasmons with wavevector k and energy 
co>k, in spite of the fact that the number density of the plasmon jV£ is explicitly absent. 

Hence it follows that for the first way only the replacement in ( |6.5| ) 

Al a (k') --^(fcV™(^U (0 U (0) ), - 4 0)c (U • • • 

gives desired contribution (similar for the conjugate term). This leads (|6.5| ) to the expres- 
sion 

-(4° )tc >3)4 o) 1^)4 o) "(^)4 o)e (^))(4° )td >04' )te >2)) (6.7) 

3 

xS(k — k\ — k 2 — k 3 )6(k' — k[ — k' 2 — k' 3 ) dkidk[ — (a <-> b, k fc', compl. conj.)}. 

8=1 

In the second case, at first step it should be replaced by 

At B (Jfe') ^ - *vl^k').P T{ - 2 ^' {A^\A^). 

This gives 

i / dk'{f^f^r v ' *v^,{ki)^f°^ 

-(A^(k[)A^(k' 2 ))(A^ 

+ (4° )t6 '(^)4' )tC '(4))K(^))(^(A; 1 )^(^))) (6.8) 

x5(k' — k[ — k' 2 )5(k — k\ — k 2 — k 3 )dkidk 2 dk 3 dk[dk' 2 + (k k' , a b, compl. conj.)}. 

By virtue of the stochasticity of gauge fields, the last term inside the parentheses vanishes. 
We replace the remaining potentials of the interacting fields by the rules 

Al{h) ^-*V uu/ (k 3 )J T ^'(A^,A^), A{{h) - Af )d (h), A e a (k 2 ) -> A^(k 2 ), etc.. 



By inspecting the expression ( |6.8j ), one sees that, complex-conjugate and non-conjugate 
amplitudes of free fields enter by a non-symmetric fashion. As will be shown below, in 



this case 5-functions (6A) type expressing the energy and momentum conservation laws 



in the plasmon-plasmon scattering have not arisen. 

By virtue of the fact that we restrict our consideration to the derivation of the kinetic 
equation for colorless excitations, i.e. 

(A® d A® e ) ~ 6 de , (A^ d 'A^ e ') ~ 6 d ' e ', (6.9) 
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all terms inside the parentheses of the expression ( |6.7| ), excepting the first term (the 
sixth-order correlator) vanish because 5 dc , 5 dc , . . . , are contracted with antisymmetric 
structure constants. Further, it is necessary to decouple the averaging of six potentials 
into pairs. We define below, what the correlations decoupling is responsible for in the 
four-plasmon decay processes. By using the definition of the correlation function ( |3.24| ). 
( ggg ), we have 

(A^\K)A^ d \k[)A^ (6.10) 

= {A^ c \k' 3 )A^ d \k[)){A^ +... 
= I u ,y(k' 3 )5 c ' d '6(k' 3 + k[)I a , x (k 2 )5 e ' d 5(k 2 - k 1 )I va {h)5 ce 5{h + k 2 ) + . . . . 

After substitution of the first term on the r.h.s. ( |6.10| ) into (|6.7| ) and performing 
integration over dkidk 3 dk' 2 dk' 3 and elementary color algebra, we obtain 

5 ab N% (~^)J { *^/t'(^)^fc^fc 2 ,-A:i^i M fc^'fe,'-^ I v i\i{—k'^I a ix{k)I va {—k 2 )dk' 1 dk 2 —{compl. conj.)}. 



As can be seen from the last expression, the required ^-functions ( |6.4| ) are not appeared 
and therefore this expression is not associated with the plasmon-plasmon scattering and 
it should be dropped. This is a general rule. The decomposition of averaging of free 
field amplitudes into the correlators containing the pair of complex-conjugate potentials 
or one of the nonconjugate potentials between the inside of the angular brackets (statis- 
tical averaging), does not give a contribution to the process of interest to us. For this 
reason, it is necessary to fully drop all contribution in the process of the plasmon-plasmon 
scattering defined by the expression ( |6.8| ), since making the correlation decoupling, the 
pair with complex-conjugate amplitudes or without conjugate necessarily arises. We write 
out decoupling of the sixth-order correlator, which gives a contribution to equation ( |6. 1| ) . 
Suppressing color and Lorentz indices and employing a condensed notion, A x = A^ki), 
we have 

(444^) = 3{(4i 3 )(4i 1 }(44) + (4a 3 ><4^2>(4a> 

+(A\,A 1 ){A\ i A 3 ){A\ i A 2 ) + (AlA 1 )(A\ l A 2 )(Ai,A 3 ) + (Al l A 2 )(A\,A 3 )(Al,A 1 ) 

+ {Al,A 2 ){A\,A 1 ){4,A 3 )}. (6.11) 

Now we consider the terms with S-functions on the r.h.s. of equation ( |3.26| ) and here, 
we also write them separately 



-i J ^'{/ bc Xtfe(4 a ( fc/ )^(^i)4(^))5(A:-A: 1 -A; 2 )dA ;i dA; 2 

- /^S^,fa<4(*)^(*0 A A - k i - k 2 )dhdk 2 }. (6.12) 
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According to the previous discussion, at a first step it is necessary to perform the 
replacement 

A^(k') — > - J tT(3)aAt '(A (0) , A (0) , A (0) ), 

A b ^(k) ^-*V^(k)J T ^'(A^,A^,A^). 

Furthermore we consequently replace the remaining two potentials of the interacting fields 
in the correlators (|612D by A°(k x ) -> - *V w ,{ki) J T ^ CV '(A^\ A®), A{(k 2 ) -> A® d (k 2 ), 
etc. This automatically leads to symmetry of contribution in A^ and A^°\ As result we 
obtain 



1 [ I *T) Jf (h'\ f bC f ff de f aC 'd fdd'e' I qlivp \y\^'u'\'<T' *jj ( t, I h. \ 

2 J t )J J J J l D fc,fc 3 ,fci + fe 2 ^fe.fcl+to^S^fc'.fe'i.fc^^ pp \ 1 ' 2/ 

-<^ 0)c '(4)4 ( < ^ (6.13) 

3 

x<5(fc — ki — k 2 — k 3 )5(k' — k[ — k' 2 — k' 3 ) dkidk' i — (a <-> 6, fc <-> fc', compl.conj.)}. 

8=1 

Because of ( |6.9| ), in the last expression it is necessary to retain only the sixth-order 
correlator. Finally adding ( |6.7D and (|6.13|) , we obtain the following equation, instead of 



a <9J af> 

—[^-(i + r^-n^fc)]^ 



: J dk'{f^f^r' 9 f 9d,e ' *Vl,(k>) fZ^M^jfi^ (6.14) 

x((41 0)c >3)4 ( ° )d >04 ( ' 0)e >D)(4 0)c (fe)4 0)d (^)4 0) lfe))^ 

3 

:5(A; — k\ — k 2 — k 3 )5(k' — k[ — k' 2 — k' 3 ) dkidk[ — (a <->• 6, A;', compl. conj.)}. 



i=i 



7 HTL-amplitudes 

Let us transform equation (6.14) into a suitable form for our further research. For this 



purpose we perform the symmetrization of the coefficient functions in the integrand on the 
r.h.s. of equation Q6.14 ) over possible permutations of color and Lorentz indices, and ar- 
guments of potentials of gauge fields within of the two groups (A^ d (ki)A^ e (k 2 )A^ c (k 3 )) 
and (Ay d (k'-^A^* (k' 2 )A„, (kQ) inside the statistical averaging angular brackets. 

For example, for the first group this symmetrization leads to the expression 

fbcf ffde ypu\a , (0)d , , x . (0)e / 7 \ a (0)c (u \ _ _ ( f bcf f fde /y^uXu yfiuaX \ 
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"rJ J \^k,k 3 ,k 2 ,ki ^k,k 2 M,kil J J K^kMMM ^k,k 3 M,k 2 ) J 

rr/yi/i^Acr -^i_Lua\ \ | (spfiaXv spfiauX \i fbcf ffde (n -r \ 

— 3! ,fel ' fc2,Jfc3 k ' k2 > kl > k z' ^ ^kMMM ^k^MM'iJ J \'- L ) 
In writing the last equality, we have used the relation between the structure constants 

yfee/ y/dc _ jbcf nfde _ jbdf jfce 

In a similar way we transform the coefficient f ac ' g f gd ' e '^yy X y y in front of the second 
group of potentials of gauge fields. 



The expression ( [71] ) is convenient because it enables us to rewrite the functions preced- 
ing the correlators on the r.h.s. of the kinetic equation ( |6.14|) in terms of HTL-amplitudes 

[9, 27, 28]. Actually, by the definition of the S-function (jl?) the coefficient in front of 
jbcfjfde equalg 



^fiuXa 
J k,ki,k 2 ,k 3 



— £ 



flUcrX 
k,k 2 ,ki,k 3 



T) ( h _l_ h \ ( Q^P Q^ pL/ \ ( caXrr qaaX 

i^paX^l -r r"2){dk,k3,ki+k 2 D k,k 1 +k 2 ,k 3 )\°k 1 +k 2 ,k 1 ,k 2 D k 1 +k 2 ,k 2 M 



+ {v <-> cr, k 2 <-> fc 3 ). (7.2) 

Furthermore, we use initial definitions ( |3.23| ), ( |3.19| ) and (|3.21| ). We present the integration 
measure d A p as dp°|p| 2 d|p|df2, where dQ is the angular measure. Using the definition of 
the equilibrium distributions ( |3.8| ) (for ji = 0) and taking into account 



+oo +oo 

N c J |p| 2 d|p| J p dp { 



ff°(Po) 
dp 



j3_ / Upl 



we perform the integral over dpo and the radial integral over d|p| in the expressions for 
S( /7 )-function ( gl9|) and S^ 7 ) -function (|3T2l| ). 



This enables us to present the expression ( |7.2j ) in the following form 

-g 2 { *r^(k, -k 2 , -h, -a*)- ^(h+ki) *r^{k, -k 3 , -h-k 2 ) *r* yff {k 1 +k 2 , -h, -k 2 ) 



- *V pa (ki + k 3 ) *T^(k, -k 2 , -k, - k 3 ) T^(fc! + k 3 , -h, -k 3 )} 



-g^T^^-k^-ku-ks), 



(7.3) 



where 



*T^ Xrj (k, h, k 2 , k 3 ) = T^ CT + ST^(k, k 3 , h, k 2 ) (7.4) 
is the effective four-gluon vertex, which represents a sum of a bare four-gluon vertex 

Y^Xa _ 2g pv g Xcr — g^g 1 " 7 — g pc ' g vX 
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and a corresponding HTL-correction 

dtt v p v v v x v a r 1 / uj 2 



5T^(k,k 1 ,k 2 ,k 3 )=m 2 D J 



47r vk + ie *-v(k + k\) + ie^vk 2 — ie 

— r~^~) - -TTTTYT-i-lT- " "T"^ - )] . V = tt. v), = 0); 

vk 3 — ie J v{k + k 3 ) + le^vki — te vk 2 — ie n 

*T^ p (k, h, k 2 ) = T^ p (k, h, k 2 ) + ST^P(k, h, k 2 ) (7.5) 
is the effective three-gluon vertex. It also represents a sum of the bare three-gluon vertex 

T^(k, k u k 2 ) = g^(k - hY + g up (ki - k 2 y + g»»(k 2 - k) v (7.6) 

and corresponding HTL-correction 

ST"(k, h,k 2 ) = ml I f ^ - , (7.7) 

j Att vk + te x vk 2 — te vk\—ie' 

m 2 D = 3Up t is the Debye screening mass. 

The polarization tensor fl3.17| ) in this notation takes the form 

dfi v»v v \ 



W v {k)=rnlU g M -u J 



Att vk + ie 



In a similar way, the coefficient in front of the product of structure constants f hd f ff ce 
in ( |7. 1| ) may be presented as 

-9 2 { *T^(k, -h, -k 3 , -k 2 )- *V pa (k 2 +k 3 ) *T^{k, -h, -k 2 -k 3 ) *T aua (k 2 +k 3 , -k 3 , -k 2 ) 

+ *V pa (k 1 + k 3 ) *T^{k, -k 2 , -h - k 3 ) *Y^{k, + k 3 , -h, -k 3 )} 

= -g 2 *t^(k,-h,-h,-k 2 ). (7.8) 

Performing transformation of the coefficient preceding (A^ d (k[)A^ e {k' 2 )A^ c {k' 3 )) 
in a similar manner, we can cast the equation ( |6.14|) in the following form 

r) dl ab 

\k 2 g^ - (1 + r 1 ) W - Tl H ^{k)}— ^ (7.9) 



• 9 
i- 



J dk' [{f bcf f fde *T^ Xu (k, -k 2 , -h, -k 3 ) + f bdf f fce *f^(k, —hi, -k 3 , -k 2 )} 



dk\ dx x 

t_ 

(3!) 2 

x *vi,(k>){r'°f° d ' e ' *t^°' x '»'(k\ -y 2 , -k' 1} -k ! 3 )+r^r' e ' *t^ x '^'(k, -k[, -4, -4)} 

3 

x5(k — k\ — k 2 — k 3 )5{k' — k[ — k' 2 — k' 3 ) dkidk\ — (a <-> b, k <-> k' , compl.conj.)]. 

i=l 
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The r.h.s. of equation ( |7.9| ) has a non-trivial color structure that actually is well repre- 
sented by non-trivial color structure of the initial dynamical equation ( |3.4|) . As will be 
shown bellow, this leads to a qualitative distinction between the elastic scattering prob- 
ability w(k, ki; k 2 , k 3 ) of colorless plasmons in a hot QCD plasma and a similar one of 
plasmons in a hot QED plasma [24, 30]. 

At the end of this Section we present the identities analogous to the effective Ward 
one in hot gauge theories [27, 28, 9]. It can be shown that the following equalities hold 

k^Y^{kMMM = *T^(k u k 2 ,k + k 3 )- *T^(k + kx,k 2 ,k 3 ), (7-10) 

k lu *T^(k, k u k 2 , k 3 ) = *T^(k + h, k 2 , k 3 ) - *F^{k, h + k 2 , k 3 ) 
(similar contractions with k 2 \, k 3a ), 

k„ *T pup {k, kx, k 2 ) = *V- lvp (-h) - *V- lvp (-k 2 ) (7.11) 

(similar contractions with k lv ,k 2p ). Here, *V~ l ^{k) = P pu (k) *A- lt (k)+Q pu (k) *A- ll (k) 
is the inverse propagator without the gauge fixing term. 



8 The kinetic equation for plasmons 

We are now in a position to explicitly compute the probability of plasmon-plasmon scat- 
tering. First of all, we make the correlation decoupling of the sixth-order correlators on 
the r.h.s. of equation ( [T.9|) in terms of the pair ones by the rule ( |6.11|) . After cumbersome 
calculations, the r.h.s. of equation (|7.9j) can be written as 



g 4 Im ( *2V(*0) / {f bcf f fde *f^(k, -k 2 , -k x , -h) + f bd ^f^ *£>*"> {h, -h, -k 3 , -k 2 )} 

x {focgfgde.fWo'Mfa _ k ^ + jadg jgce *f t^'AW^^ _ k ^ _ k ^ _ k ^ (g ^ 

xlw'{h)hx'{h)laa'{k 2 )8{k -kx-k 2 - k^dk^^d^. 
In deriving of ( |8.1| ) we have used two relations, which satisfy *T puXcr (k, k\,k 2 , k 3 ) 
*f p ^(k, h, k 2 , k 3 ) + *f pX ^(k, k 2 , kx, k 3 ) + *f^\k, kx, k 3 , k 2 ) = 0, 

*f^(k, kx, k 2 , k 3 ) = *f fiaXv (k, k 3 , k 2 , kx). (8.2) 

Their correctness may be verified by a direct calculation using known properties of HTL- 
amplitudes [27, 28], entering into the definition of *T puXa (k, kx, k 2 , k 3 ) ( [7T3p . The first 
of the relations in (|8.2|) supplements, in the some sense, a similar relation for the HTL- 
correction dT^ to the bare four-gluon vertex, originally proposed by Frenkel and Taylor 
[28]. The second one is a property of invariance of the function *T, when the momenta 
order is reversed. Notice that order of the space-time indices in equation (|8.2| ) is important. 
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Furthermore, if we restrict our consideration to the study of plasmons scattering by 
plasmons, then in the spectral decomposition of I^i into the orthonormal projectors 
and Qnn', it is necessary to keep only the longitudinal part, i.e. set 

hw(h) = Q^(k 3 )I l (h) , Jaa'(*i) = Qx\>(ki)I l (h) , I aa ,(k 2 ) = Q aa ,(k 2 )I l (k 2 ). (8.3) 

In the propagator *V^(k) we also keep only the longitudinal part — *A\k). To 

determine the imaginary part of *A (k) = l/(k 2 — n'(A;)) we use the approximation (see 
e.g., Pustovalov and Silin in [22]) 

in sign (Im (k 2 - Yi\k))5{Re (k 2 - Tl\k))) = 

signu5{Ree l {k)), (8.4) 



k 2 -W(k) Re(k 2 -W(k)) 

P in 



k 2 Ree l (k) k 2 

which holds for a small lme l (k). Here, we consider that because of ( |2.3| ), sign (Ime l (k)) = 
signu. The symbol P denotes a principal value. We perceive the 5-function of the real 
part of the longitudinal permeability, which appears in ( |S.4| ) in the ordinary sense ( |6 . 6| ) . 
Substituting ( |6.6|) into ( |3.4|) and choosing uj = uj 1 ^ > for defmiteness, we obtain the 
required relation 



Im*2V(*0 « ^^{ ^li^ y 1 .^-wL). (8-5) 



With allowing ( |3.3|) , ( |8.5|) and performing the color algebra with the following identities 



jbcf jfde jacg jgde _ pj2 ^ 



ab 

C w 5 



jbcf jfde jadg jgce _ _ jyS^ab 



we rewrite ( |8.1| ) in the form 
1 r dRee l (k)^-i 



|2 



+ Re ( f (*, -fa, -fa, -fa) f •(*, -fa, -fa, -fa))} mk)mki ^ (k2mk3) (8-7) 
xl l (k 1 )l l (k 2 )l\k 3 )6(k -k x -k 2 - k z )dk x dk 2 dk z . 

Here, we denote 

*f(k, h, k 2 , h) = *f^(k, k 1} k 2 , h)M k )Mh)uxfa)uM. (8-8) 

By virtue of the fact, that the spectral intensity functions I l (ki),i = 1, 2, 3, entering into 
satisfy the equations Kee (ki)I (ki) = [10], they have the following structure 

I 1 (fa) = I&fa - O + IV(<* + O , z = 1, 2, 3, (8.9) 
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where !{. are certain functions of the wavevectors kj. 



We substitute ( |8.9|) into (|8.7| ) and perform the integration over the frequency dcuj with 
the help of the ^-functions. The function I h in the l.h.s. of the kinetic equation ( |7.9| ) also 



has a structure of form (|8.9|) . Here, we keep only a positive branch u > in agreement 
with flS.5| ). Furthermore performing integration ( |7.9j ) over du, where the r.h.s. of ( |7.9| ) 
has the form (|8.7|) , as result we obtain 

u2f dRee l (k) ^ dll l2 ^ dRee l (k) ^ dll 

= ^ 4 ^( d R ^ (fc) )~I 4 / fci, **, W* + k 1 + k 2 + ha) 
+Q(k, -ki, -k 2 , —k 3 )8(k - k\ - k 2 - h) + Q(k, k 1} k 2 , —k 3 )5(k + h + k 2 - k 3 ) 

+Q(k, fci, -k 2 , k 3 )5(k + k x -k 2 + k 3 ) + — fci, fc 2 , fc 3 )5(A; - fc x + fc 2 + k 3 ) 
+Q(k, fci, — A; 2 , -k 3 )5(k + k\ — k 2 — k 3 ) + — fci, fc 2 , -k 3 )5(k — k\ + k 2 — k 3 ) 
+Q(k, -h, -k 2 , k 3 )5(k -k 1 -k 2 + k 3 )} on ^ she iill x lljl 3 dki dk 2 dk 3 , 
where we have defined 

W ' **' fe) S u^u^u^u^h) (8 ' U) 

x{| *f (*, fc 2 , k u k 3 )\ 2 + | *f(/c, h, k 3 , k 2 )\ 2 + Re ( *f(k, k 2 , k u k 3 ) *f fcx, fc 3 , fe))}. 

In equation ( |8.10| ) the z-dependence of ii, etc is understood, although not written 
explicitly. 

The first term inside the square brackets on the r.h.s. of flS.lOp describes the process 
of simultaneous fusion or emission of four plasmons in the plasma. Considering that 
the 5-function has no support on the plasmon mass shell, its contribution to the kinetic 
equation vanishes. The second, third and fourth terms describe the decay of one plasmon 
into three or the fusion of three plasmons into one. As was mentioned in section 2, these 
processes are forbidden by the conservation law ( |2.23| ). Therefore these terms also vanish. 
The remaining terms describe the scattering of two plasmons off two plasmons. Three 
terms imply the existence of three possible channels of a given process, which change the 
plasmon numbers density iV£ 

g* + g*l ^ g2 + §3> g* + g2 ^ gl + g3> g* + g3 ^ Si + g2- 



If we perform replacements ki «-> k 2 in the next to last term on the r.h.s. of (|8.10|) , and 
kx <-> k 3 in the last term, then the r.h.s. of fl8.1U|) can be presented as 



ng*N 2 c ^( dRe d £ ^ k) )^ J [Q(k, k u -k 2 , -k 3 )+Q(k, -k 2 , k u -k 3 )+Q(k, -k 3 , -k 2 , h)} 



on— shell 
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xS{J k 



4 2 - 4 3 )5(k + kx - k 2 - k 3 )44 2 /[ 3 dk 1 dk 2 dk 3 . 



By using the definition of kernel Q ( |8.11| ) and the properties ( fj.2p , it can be shown, that 
Q(/c, -& 2 , fa, -fa) + Q{k, -fa, -fa, fa) = 2Q(k, fa, -fa, -fa). 



Taking into account the last relation, it goes over from the function to the number 
density of plasmons 



Ni 



,dRee l {k) 



I 1 



and recovering the complete form for finite values iV k , we obtain the required Boltzmann 
equation ( |2.25| ), where the probability of plasmon-plasmon scattering is defined by the 
following expression 



W (k,k i; k 2 ,k 3 ) = 37r(27r)ViV ( 2 



x 



sdRee l {fa)^ 


-i/dRee l (fa)\ 


^ duj2 ' 





-1 



+Re ( *f (fc, -fa, fa, -fa)*P(k, fa, -fa, -fa))}- 



1 , dRee l (k) yi, dRee l (fa) 

k 2 k\k\ki v ' \ 

{| *f (/c, -fa, fa, -fa)\ 2 + I *T(k, fa, -fa, 

1 



(8.12) 



-fa)\ 2 



v? ( k) v? ( fa ) v? ( fa ) u 2 ( fa ) 



on— shell 



The result ( 8.12|) is rather unexpected. As we can see from the expression ( 8.12|) , this 
probability does not reduced to the squared module of one scalar function, as this occurs 
in the Abelian plasma. Here, the scattering probability is defined by the squared module 
of two independent scalar functions and their interference?], in spite of the fact that in this 
paper we restrict our consideration to a study of the nonlinear interaction of only colorless 
excitations. This radically distinguishes the Boltzmann equation (|2.25|) , describing the 
effects of the collisions among colorless soft excitations from the corresponding Boltzmann 
equation including the effects of the collisions among colorless hard excitations [34, 18, 19]. 
In the last case, the Boltzmann equation, corrected to color factors, fully coincides with 
the corresponding one in the Abelian plasma. This is a point which we find difficult to 
interpret and therefore additional analysis of this problem is required (see, also discussion 
in conclusion). 

The scattering probability can be written in the form, which is manifestly symmetric 
under permutations of the external momenta k, k 1; k 2 and k 3 



^(k,k i; k 2 ,k 3 )=7r(27r)ViV ( 2 



. rv rv^ I\j2 



s dRee l {k) yis dRee l {fa) yisdRee l {fa)^-i 



duo 



x 



'<9Ree'(A;3 



{\ *r(k, -fa, fa, -fa)\ 2 + | *f {k, fa, -fa, -fa)\ 2 + 1 *t(k, fa, -fa, -k 3 )\ 2 + 



It is easy to see from the expression (8.12), that the function w(k, ki; k2, k/j) is positively definite. 
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Re ( *f(k, -k 2 , k h -k 3 ) *f+(/c, h, -k 3 , -k 2 ))+Re ( *f(k, -k 2 , h, -k 3 ) *f t(fc, h, -k 2 , -k 3 )) 

1 



+Re ( *T(k, h, -k 3 , -k 2 ) T\k, h, -k 2 , -k 3 ))}- 



on— shell 



u 2 (k)u 2 (ki)u 2 (k 2 )u 2 (k 3 ) 

This expression is suitable for checking the symmetry conditions (|2.26|) , which is imposed 
on the plasmon-plasmon scattering probability. 

Recall that the function T, which appears in the expression for probability ( |3.12|) , 
is defined by expressions ( fc>.8| ) and (|7.3|). As we have shown in [10], the expression ( |S.8| ) 
(exactly, its part, independent on a gauge-parameter) is gauge invariant at least in a 
class of covariant and temporal gauges. With regard to the reasoning in section 4 this 
automatically leads to gauge invariance of the kinetic equation ( |2.25| ). The problem of the 
dependence of the plasmon-plasmon scattering probability ( |S.12D on a gauge parameter 
coming from the gauge fixing term in the gluon propagator, is discussed below. 

Let us separate the terms with a gauge parameter from the expression *T(k, —k 2 , k\, —k 3 ). 
By using the Ward identities ( 7.10 ) and 1 1| ) , we have 



t{D P a{-ki + fc 2 )A°(-A; 1 + k 2 ) *T^(k, -k 3 , h - k 2 ) T^-Aa + k 2 , k u -k 2 ) 

+D pa (-h + fc 3 )A°(-A:i + k 3 ) *r^(k, -k 2 , h - k 3 ) T^-fci + k 3 , k h -k 3 )} 

x u^(k)ux{ki)u a (k 2 )u v (k 3 ) (8.13) 

= £{(A°(-^i + k 2)f{ *V^^(k) - *V~ llMJ {k 3 )){ *V~ 1Xv {k 2 ) - *V- 1Xv (-h)) 

+ (A°(-fci + k 3 )) 2 (*V~ l ^(k) - *©- 1 ^(A; 2 ))(*I?- 1Aff (A;3) - *V~ l Xa (-h))} 

xU p (k)ux(h)u a (k 2 )u u (k 3 ). 

It is easily shown that expression on the r.h.s. of ( |3.13|) vanish either because *V~ ltiiy (k) 
is transverse, or by the definition of the mass-shall condition, i.e. 

k^V-^ik) = 0, *V- 1Q \k)l =< = 0. (8.14) 

If we carry out the following replacements on the r.h.s. of the expressions (^.8|), ( |7.3| ) 

k 2 

u^k) -> u p (k) = j^(ku - u u (ku)), 



and 



"DM - *v pa {k) = *v pa {k) - (^j^c pa (k) + ^Ldm) * Al (*) 



-iD pa (k)/\\k)-i Q7 —D pa (k), 



A 2 

' (ku)' 

where £ is a g au g e parameter in the temporal gauge, then it can be proved that a similar 
statement holds in the temporal gauge also. 
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Thus, the gauge-dependent parts of w(k, ki; k 2 , k 3 ) drop out, since they are multiplied 
by the mass-shell factor. These factors are proportional to {u — wy. However, in the 
quantum case Baier et al [35] observed that a naive calculation of the gluon damping 
rate in a covariant gauge appears to violate this consideration. The mass-shell factor is 
multiplied by the integral involving a power infrared divergence which is cut-off exactly 
on the scale (a; — ~ g 2 T. This problem was careful discussed in consideration of the 
nonlinear Landau damping [10], when in (|7.3|) and consequently ( E>.13|) it is necessary to 



set 

ki = —k, ki = —hi, k 3 = k\. (8.15) 

We have shown, that the gauge-dependent part of the nonlinear Landau damping rate 
vanishes on a mass-shell at least for the plasmons with zero momentum. In the case 
of the process of the plasmon-plasmon scattering considered here, the problem on the 
dependence of w(k, k^ k 2 , k 3 ) on a gauge parameter is more subtle, since instead of flS.15| ), 
we have condition: k+k± = k^+k^, only. Research into this nontrivial question goes beyond 
our present goal. 



9 Lifetimes of plasmons 

To calculate the lifetimes of colorless plasmons we first linearize the Boltzmann equation 
( p. 25| ) (here, we assume that off-equilibrium fluctuation is perturbative small), writing 
the number density of plasmons as 

Nl = N l eq (k) + 5Nl 

where iV^k) = (e w k/ T * — 1) _1 is the Planck distribution function and T* is a certain con- 
stant, which can be interpreted as a plasmon gas temperature in the statistical equilibrium 
state. Then we find 

= / (^ (^ ( 27r ) 4 ^ + < - < - <W k + kx - k 2 - ksMk, k i; k 2 , ka) 
x{5iV[[iV^(k 2 )iV^(k3)«(k 1 ) + 1) - ^(k x )(iVi g (k 2 ) + l)«(k 3 ) + 1)]+ 
5Nl 2 [Ni q (k 3 )(N l eq (k)+l)(N^ <-> kx, k 2 <-> k 3 )}. 

This equation can be further simplified if we use the following parametrization for off- 
equilibrium fluctuation of the occupation number 5iV k [12, 18] 

diVyk). 

i 



SNl = = (l/T*)N l eq (k)(N l eq (k) + l)Wl (9.2) 
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Let us introduce the momentum transfers, which carries the soft gluon exchanged in the 
collision of two plasmons, setting 



k 2 = k - q, k 3 = kx + q. (9.3) 

Performing the integration over dk 3 , replacing the integration over dk 2 by one with re- 
spect to momentum transfer and taking into account ( |9.2| ) and ( |9.3| ), finally we derive a 
linearized Boltzmann equation for W k function 

dWj , y , dWj f did f dq JV^KiV^k-qj + lKiV^ kx + ^ + l 



dt k <9x " J (2vr) 3 J (2vr) 3 (A^(k) + 1) 

x u>(k, k i; q) (2n)5U - 4_ q + < - < +q ){W k - W k _ q + W kl - W kl+q }. (9.4) 
Here, we goes over from the function u>(k, ki; k 2 , k 3 ) (|8.12|) to a new function u>(k, ki; q) 

u;(k,ki;k2,k 3 )|k 2 =k-q,k 3 =k 1 +q = u>(k, ki;q). 



Based on the exact form of the r.h.s. of equation (|9.4| ), we define the lifetimes of the 
plasmon of momentum k as follows 

1 r dk a f dq NiMiN^-^ + l^NUk. + cO + l) 



r pl (k) J (2vr)3y (2tt)3 W + l) 

x w(k, ki; q) (2tt)(5(4 - 4_ q + J kl - 4 1+q ). (9.5) 

Here, the integrand has a more involved structure in comparison with a similar expression 
for the case of the lifetimes of the hard transverse gluon [34] . The reason of this fact is 
that the momentum of the soft quasiparticles (plasmons) becomes of the same order as 
momentum transfers and under these conditions one must take into account the non-trivial 
character of the frequency dependence uj of the collective gluon modes on momentum k 
and vertex corrections. 



It is not difficult to estimate the order of the expression ( p.5Q at the momentum scale 
gT. Considering the plasmon gas in thermal equilibrium with hard particles from the 
heat bath, i.e. T* ~ T, and using the definition of w function, we obtain 



— ~ g z N 2 c T. (9.6) 

Tpl 



However, obtaining a numerical factor of proportionality in ( |9.6| ) is a complicated problem 
even for limiting case of k = 0-mode. Here, we restrict our consideration to following 
general remark, which somewhat simplifies the matter. 

In proving the gauge invariance of the nonlinear Landau damping rate, we have shown 
[10] that the *T function (|8.8|), ( [7.31 ) entering into the definition of the scattering proba- 
bility (|8.12 ) can be introduced in its simplest form 

*f(k,-k 2 ,k 1 ,-k 3 ) = 
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k 2 k\k\k\{ *T 0000 (k, -k 2 , h, -k 3 )- *V pa (k-h) *T 00 ^k, -k 3 , h-k 2 ) *T a00 (-h+k 2 , h, -k 2 ) 

- *V pa {k - k 2 ) *T °P(k, -k 2 , k, - h) *T am (-ki + k 3 , h, -h)}\ on _ shell . 

However, even if the last expression is accounted for, after performing the integration over 
solid angles in <5r 0000 and 5T 00p HTL-amplitudes (see, e.g., [28]), we obtain expressions 
which are very cumbersome. For this reason passage to the limit |k| — > is non-trivial. 

The estimate ( |9.6|) shows, that at the momentum scale gT a " collision" damping rate 
of the soft longitudinal mode is suppressed by a power of g relative to a value of the 
nonlinear Landau damping rate (~ g 2 T) [10]. Therefore at the soft momentum scale it 
can be neglected by the influence of the plasmon interactions among themselves on the 
relaxation dynamics of a plasma excitations. 

As was mentioned in Sec. 2, the process of the nonlinear Landau damping leads to an 
effective pumping of energy across the spectrum towards small wavenumbers. By the effect 
of pumping, all plasmons will tend to concentrate near small |k| = |k | — > 0. However, 
phase space, within which the plasmons are occupied, proportional to |ko| 3 , will also be 
very small (a similar state, when a great many plasmons with small wavenumbers are 
mainly concentrated in a small phase volume, in the theory of strong turbulent electron- 
ion plasma this is sometimes called a plasmon condensate [36]). By virtue of this fact 
intensive collisions between plasmons arise and this scattering process is described by 
the Boltzmann equation (|2.25| ). This leads to a "throw out" of plasmons from region of 
small |k| and thus a suppression of the increase of the k = 0-mode. In mathematical 
language this denotes that for definite values of the momentum, the magnitude of the 
off-equilibrium fluctuation <5iV£ becomes as large as N l eq (k), and therefore a linearization 
of the Boltzmann equation ( |2.25| ) is no longer valid, like the estimation ( |9.6j ), corrected 
within the framework of this approximation. In the region |k| <C gT, where collisions 
among plasmons start to play a role, it is necessary to solve the exact nonlinear integro- 
differential equation, whose r.h.s. is considered as both the process of the scattering 
plasmons off thermal particles, and the process of four-plasmon decays, i.e. 

+ vL !r = -^hf J ^ (4 - O Q(k, ko^< (9.7) 

+ / (2^ Wf Wf {2nr6U + 4l ~ 4a " 4j5(k + kl " k2 " k3) 

x W (k, k i; k 2 , k 3 ) + ;vXX - ivXX - *k«), 

where 

Q(k , kl) = (W^r 1 f^my 

x J ^(^-<-v(k-k 1 ))|^ c (k,k 1 )+^ll(k,k 1 )+>l ± (k,k 1 )| 2 . 
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The amplitudes Ai c , Ai^ and Ai 1 - have the following forms [10] 

g 2 1 (kv)(k lV ) 



^ c (k,k 1 ) = 



|k||ki| uLu^ loI- (kv) 
<? 2 ((k-kx)v) 



x 



|k||k!| (k-k x ) 2 

f ( fc - fc l) 2 ' A '( fc - fc l) JT^ -Jfc + h)k%{k - h) 1 ) 

g 2 (Kk-kOH 



X 



M ™ |k||k,| „>(k- kl )» 

*A t (k-k 1 ) 



UJUJ\ 



*T ijl (k, -k u -k + fci)FA;i[n(k - k x )] 



Here, n = [kki]. We present the amplitudes Ai c , tM" and Ai 1 - in their simplest form, 
defined by temporal gauge. 

A qualitative analysis of the solution behavior of a similar equation in the case of 
electron- ion plasma with a rather crude approximation of the kernel Q(k, ki) and the 
probability of plasmon-plasmon scattering w(k, ki; k2, I&3), was performed by Kovrizh- 
nykh [30]. He has shown, in particular, that really in the kinematic regime of small 
wavenumbers |k|, the process of plasmon-plasmon scattering starts to play a dominant 
role, that prevents an increase of the k = 0-mode. In the case of the the non-Abelian 



plasma, an analogous qualitative investigation of the solution of equation (|9.7| ) presents 
difficulties because of the above-mentioned complexity of the integrands on the r.h.s. of 
equation ( |9.7|) , and therefore it requires some additional approximations, allowing these 
functions to be made more visible and suitable for numerical calculations. 



10 Conclusion 

Within the framework of the semiclassical kinetic theory of a hot gluon plasma we have 
obtained the transport equation for the plasmons, taking into account four-plasmon de- 
cay. The probability of plasmon-plasmon scattering at the leading order in the coupling 
constant is derived. This is defined with the help of three-gluon and four-gluon effective 
vertices, and the effective propagator only, as the probability of the process of scattering 
of the plasmon by hard QGP particles [10]. It is proved that the scattering probability is 
gauge-independent at least within a class of covariant and temporal gauges. 

In this paper we have restricted our consideration to a derivation of the Boltzmann 
equation for the simplest collective Bose-modes of hot gluon plasma, colorless plasmons. 
This fact was reflected in deciding on a trivial color structure of the plasmon occupation 
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number: N{? b ~ 5 ab . It is clear that such a choice of trivial color structure is completely 
neglected by the fundamental property of plasmons in the non-Abelian plasma, such as 
the availability of color charge. This is a very important difference from the Abelian case, 
where the plasmons do not carry electric charge. While the initial state of the problem 
can have a trivial color structure, the dynamics of these collective modes can change the 
color of the plasmons. In light of these remarks, here we are dealing with an Abelianized 
version of more complicated colored plasmon dynamics. It is a first step towards a full 
description of plasmon transport properties to derive the Boltzmann equation which takes 
into account the precession of the color charge of the plasmon. It is probably only in the 
construction of such a complete equation that we gain a reasonable explanation of the 
structure of the scattering probability (|8.12|) . 

The scheme of the derivation of the Boltzmann equation for colorless plasmons, pro- 
posed in this paper, admits the straightforward generalization to the case, when the 
number density of the plasmon iV£ af) has a non-trivial color structure, e.g. such as 
jyiab _ (y c ) afe jv£ c . i n this case, as was discussed in Sec. 4, it is necessary to assume 
that the regular (background) part of the field A^(X) is different from zero. The density 
matrix, effective propagator, three-gluon, four-gluon vertices become the functionals of 
the background field Aj^(X) [9, 18], which assumed vanishing at Xq — > — oo. This signif- 
icantly complicates the problem of construction of the required kinetic equation. It can 
somewhat simplify the problem if it is considered that the fluctuation 5Nl_ ab is a small 
perturbation in relation to the equilibrium value 5 ab N l eq (\i), where N l eq (k) is the Planck 
distribution, and restrict the consideration to the linearized Boltzmann equation for col- 
ored plasmons. The exact consideration of a given question will be subject of separate 
research. 

There is an independent test of the validity of the derived Boltzmann equation ( |2.25| ) 



(exactly, its linearized version ( |9.4|) ). As was mentioned in introduction in [18] the funda- 
mental derivation of the Boltzmann equation for a high temperature Yang-Mills plasma 
was proposed by Blaizot and Iancu within the CTP formalism framework. Their deriva- 
tion relies on a gauge- covariant gradient expansion of the Kadanoff-Baym equations for 
the gluon two-point function. The Boltzmann equation has emerged as the quantum 
transport equation at leading order in g for the gauge-covariant fluctuation 5G of a hard 
gluon propagator. 

Besides the above-mentioned paper, the Kadanoff-Baym equations for the off-equilibrium 
propagator of the soft gluon V^X^Y), which are formally identical to those for a hard 
gluon propagator G^ U (X, Y), are written out. These equations are used in [18] only 
to deduce the relation between the off-equilibrium gauge-covariant fluctuation V < (k,x) 
and the gauge-covariant fluctuation of the leading-order soft gluon polarization tensor 
511 < (k,x), and the problem of self-interactions of the soft fields is not considered here. 
However, in principal there is nothing to forbid the use of these equations for research in 
the dynamics of the soft fields and construction of the relevant transport equation within 
the framework of the scheme, suggested by Blaisot and Iancu [18]. Here, by £]!<(£;, x) 
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we mean the fluctuation of the next-to-leading order of the soft gluon self-energy, involv- 
ing three- and four-gluon off-equilibrium vertices with soft external lines. The relevant 
effects of self-interaction of the soft fields are encoded in these functions. Here, we note 
that in the real-time formalism, technical complications resulting from the doubling of 
degrees of freedom are arisen. However in recent years a suitable technique allowing one 
effectively to work not only with the single-particle propagator, but also with three- and 
four-point functions, has been developed [8]. This greatly simplifies calculations in real 
time and enables us to derive the transport equation for the soft gluons, in particular, 
the plasmons, directly from the underlying quantum field theory and compare it with the 
equation, obtained in this paper in the context of the semiclassical approximation. This 
formal scheme is also needed to specify the limits of validity of the semiclassical kinetic 
approach to the research of the processes of nonlinear interaction of the soft fields in a 
hot QCD plasma. 
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